Exhibition of the periodicity of Quantum Fourier Transformation 
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^ ■ Abstract 

O 

(N 

The remarkable capability of quantum Fourier transformation (QFT) to ex- 

^. 

Q_i, tract the periodicity of a given periodic function has been exhibited by using 

+-> ' 

nuclear magnetic resonance (NMR) techniques. Two separate sets of exper- 
ts 

2 . 

q-( iments were performed. In a full QFT, the periodicity were validated with 



state tomography and fidelity measurements. For a simplified QFT, the three- 



qubit pseudo-pure state was created by introducting an additional observer 
spin, and the spectra recorded on the observer spin showed intuitively the 
power of QFT to find the periodicity. Experimentally realizing the QFT pro- 
vides a critical step to implement the renowned Shor's quantum factoring 
algorithm and many other algorithms. Moveover, it can be applied to the 
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study of quantum chaos and other quantum information processing. 
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The most striking discovery in quantum computation is that quantum computers can 
efficiently perform some tasks which are not feasible on a classical computer. For example, 
the most famous of the current quantum algorithms — Shor's quantum factoring algorithm 
|jlj], can factor a number exponentially faster than the best known classical algorithms, which 
leads to the cracking the RSA encryption system |2|. However, the key to Shor's algorithm 
rests in an essential way on the application of a suitable quantum Fourier transform to 
efficiently determine the periodicity of a given periodic function ||. The QFT is the key 
ingredient for not only Shor's algorithm but also many other interesting quantum algorithms, 
including Simon's algorithm (Simon 1994) fU|, Kitaev's algorithm |||7|, the estimating 
arbitrary phase algorithm |||8| and the ordering-fmding algorithm || etc.. Moveover, the 
scheme of exploiting quantum computation to study quantum chaos proposed by Schack || , 
Brun and Schack |1(J, are also based on the QFT. 



For an arbitrary positive integer x, the QFT can be defined by a unitary transformation 

i N-l 

QFT : \x) - -= £ e 2 ™v/ N \y), for x G {0, 1, 2, T - 1} , (1) 

where N = 2 n , with n being the qubit number. Recently the QFT have been realized 
on a three bit NMR quantum computer. However, the implementation of the QFT was 
performed on the thermal equilibrium state |13| so that the striking power of the extraction 
of a periodicity wasn't exhibited from the experimental results. 

Imagine an oracle that computes a function / : Z n — > Z m that has an unknown period r 
(1 < r < 2 n ): 

f(x + mr) = f(x), (2) 

where m is any integer such that x and x + mr lie in {0, 1, 2, 2 n — 1} . The goal is to find 
the period r of f(x). Classically, this problem is hard, requiring O (2 n ) queries to hit two 
equal values with high probability [jTJJ]. But a quantum algorithm described as follows can 
find r in time poly(n) [|12|.|15|]. 



Define a unitary transform U : 

U\x)\y)^\x)\y®f{x)), (3) 

where © denotes addition modulo 2. One first initialize two registers to |0) |0) , then create 
an equal superposition -^=J2 X \ x ) i n the first register (easily prepared by applying H^). 
Owing to quantum parallelism, querying the oracle and applying U prepare the state 

1 N-l 

l/) = 7^EWI/W>' ( 4 ) 

By measuring the second register, the result \f (xo)) f° r some < xq < r — 1 is obtained .In 
the meantime, the first register is collapsed into the state 

I K-l 

J2 \ x o+jr) , (5) 



K j=0 



where N — r < xq + {K — l)r < N. The periodicity r is not explicitly exhibited due to the 

Offset Xq. 

Applying QFT to the state , one gets 



1 i ■w^Qj .N\ 



F\^) = ^E/ m ^ (6) 

A value c is supposed to be obtained by measuring the register, where c is necessarily a 
multiple of N/r, i.e., c/N = A/r, (0 < A < r). If A is fortuitously coprime to r, r can 
be determined by cancelling c/N down to an irreducible fraction. According to the prime 
number theorem [|T^ , p~6|JT~7[] , the probability that A is coprime to r is at least 1/logr which 



exceeds 1/logN. Hence repeating the above procedure 0(logN) times one can succeed in 
determining r with any prescribed probability 1 — e as close to 1 as desired. 

Another way of analyzing the distribution of outcomes obtained is to calculate the re- 
duced density matrix for the first register and calculate the measurement statistics. As the 
NMR quantum computer is ensemble quantum computation, we can achieve the statistical 
results of the probability distribution of the final state. It can be seen from Eq. (5) and 
(6) that, the QFT transforms the input state with the period r into the output state with 



the period k = N/r containing r items due to the displacement invariance of the QFT. For 
example, if the 3-qubit state in Eq. (5) input to QFT is of the form 



~(|0) + |2} + |4> + |6» 



or 



1 



(|1> + |3> + |5> + |7», 



(7) 



(8) 



the state periodicity should be r = 2. Applying the QFT, one gets 

1 



QFT 



V2 



(|0) + |4» 



(9) 



or 



QFT: 



V2 



(|0>-|4». 



(10) 



From the state periodicity after the QFT k = 4 one can infer the state periodicity before the 
QFT r = N/k = 2, which is consistent with the state Eq. (7) or (8). Contrarily, if has 
the form of Eq. (9) or (10), the state periodicity after the QFT k = 2 so that r = N/k = 4. 

We have demonstrated these effects by liquid-state NMR using the carbon- 13 labeled 
alanine NH£ -C a H(C l3 H 2 )~C ' dissolved in D 2 0. All NMR experiments were performed 
on a Bruker ARX500 spectrometer with respect to transmitter frequencies of 500.13 MHz 
i}H) and 125.77 MHz ( 13 C). The measured NMR parameters are listed in Table 1. We can 
see from above that, the second register is only used to prepare the periodic state in the 
first register. Thus, in actual experiments, we just used the first register for determining the 
periodicity. In addition, owing to the nature of NMR ensemble, the display of the final state 
were performed by two different methods, tomography |22j and "spectral implementation" 
23|. 



The logic network for the QFT is shown in Fig. 1 ^,11,18 , consisting of a one-qubit 



gate 



H 



3 y/2 



(ll \ 



v 1 - 1 / 



(11) 



acting on qubit j and a controlled- Rd gate 

d ^ 



R 



d 



y e-/ 2d j 



(12) 



acting on qubit j conditional on qubit k being in the |1) state with d = j — k. Here the 
Hj gate can be realized by the NMR pulse sequence Xj(ix)Yj (|) or Yj ( — f) and 
the controlled-i?j gate, by the pulse sequence Zj(^)Z k Jj k (— f ) p^9| , |20| . By carefully 
choosing the expressions for the Hadamard and Z gates and exploiting the commutation 
relations |21[| to eliminate as many unnecessary operations as possible, the complete pulse 
for the QFT can be reduced as: 



(applied from the left to the right). Finally a swap operation Sis are used to reverse the 
order of qubits to obtain the desired output from the QFT. We could just as easily relabel 
qubits 1, 2 and 3 in place of the Swap gate. Thus, we need not perform an actual physical 
swap gate here; a mental relabeling of the qubits is sufficient. 

We preformed two separate sets of experiments. In the first set, the full QFT was 
executed on a three-qubit quantum computer. Three 13 C nuclei were chosen as three qubits 
and protons were decoupled during the whole experiments by using a standard heteronuclear 
decoupling technique. All 13 C nuclei were taken 0.7ms for selective | pulses of a Gaussian 
shape. At the beginning of the experiment, we labeled C, C a and C 13 as spin 1, 2 and 3, 
respectively. After the QFT, we identify spin 1 with C 13 spin 3 with C . 

The pseudo-pure state was prepared from the thermal equilibrium state by the procedure 
summarized in Table 2, in which the magnetic field gradients (denoted by G z ) to dephase 
off-diagonal elements of the density matrix at strategic points along the way were used pi] . 



Owing to the difference of the relaxation times between nuclei, we used a ir rf pulse and a spell 
(the reversal recovery) to balance the effect of the different relaxation. Meanwhile step (1) in 
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Table 2 was replaced by the spell. By adjusting carefully the spell, the experimental spectra 
for a pseudo-pure state p oo — ^1*^2^3 (If — |0) (0| = \ + Ii z ) were finally recorded (shown in 
Fig. 2b) through reading-out pulses on the spectrometer. The normalized deviation density 
matrix p oo was confirmed by full tomography |EJ shown in Fig. 3a. 

The state in Eq. (7) and Eq. (9) can be prepared by the Hadamard gates H1H2 and 
Hi from the pseudo-pure state, respectively. Applying the sequence (13), the experimental 
spectra were shown in Fig. 2 (c) and (d). The density matrices were reconstructed from 
the obtained experimental spectra. The real parts of these matrices are shown in Fig. 3 (b) 
and (c), respectively. The fidelity of the QFT were calculated using the measure, called the 
attenuated correlation ||2"4} 



Tr(p th p exp ) 
Tri^p^p 1 



Here, p th is defined as a theoretical output state, transformed by the ideal transformation 
on a computer for the measured pseudo -pure state Pooo' &frd an experimentally implemented 
control sequence for the same transformation on the spectrometer to get p exp . The values 
of the correlation for each of the three tomographic readouts were c(pqqq) = 1 (by defini- 
tion), c(p^ 1 ) = 0.92 and c(p^ 2 ) = 0.84. The measure reflects the imperfections of the 
experiments, including inhomogeneity of RF fields and static magnetic fields and imperfect 
calibration of rotations, gradient pulses and relaxation. 

In the second set of experiments, we prepared a pseudo-pure state with introduction of 
an observer spin and employed a further simplification of the QFT proposed by Preskill 



|25fl, in which no two-qubit gates are needed at all and only n Hadamard gates and n — 1 
single-qubit rotations are applied. Its key thoughts is the symmetry of the controlled-/^ 
gate on the two qubits, and to measure a single qubit first and then apply the controlled- Rj 
gate to the next qubit, conditioned on the outcome of the measurement of the qubit, instead 
of applying controlled-/^ and then measuring. Of course, it is only suitable for a special 
state, for example, the computational basis state IX3) \%2) \x\) , (xj = or 1) and the state 
in Eq. (7) and (9). 



The pseudo-pure state IqzI^I^I^ with introduction of an observer spin 1$ were prepared 



by a spatial averaging method proposed by Sakaguchi et al. P6[ , the procedure shown in 
Table 3. We chose C a as the obssever spin Jo due to resolved scalar J couplings to all other 
spins and C, C 13 and H being joined directly with C a as Spin 1 (ii), 2 (J 2 ) and 3 (I3), 
respectively. The methylic hydrogen nuclei were decoupled using the continuous wave (CW) 
mode. 

The results of the second set of experiment were shown in Fig. 4. The transitions that 
connect the Iq and Jq manifolds of the observer spin Iq are assigned to the computational 
basis states of the other spins. The readout on the observer spin Iq can tell us which states 
the computational spins I\, I2 and 1% were in |p3[| . As exhibited in Fig. 4b, only one line 
corresponding to the 1^1% 1% state indicated the computational spins I\, I2 and J3 was in 
the 1 000) state. It can been seen from Fig. 4c that, after swapping Spin 1 and Spin 3, I\, 
I 2 and I3 was in |000) + 1 100) = |0) + |4) , thus the state periodicity k = 4 and the state 
periodicity before the QFT r = N/k = 2 was deduced; From Fig. 4d, I±, I2, I3 was in 
1 000) + 1 10) + 1 100) + 1 1 10) , hence, k = 2 and deduced r = 4. This simple readout method 
can make one quicker and more intuitive to get the results and is suitable for any number 
of qubits. Compared with the tomography JEJ used in the first set of experiments, which is 
practically feasible for no more than three qubits, it is simpler and more efficient. However, 
its signal-to-noise ratio (SNR) is not very good and it is more demanding for the sample. 

In conclusion, using the NMR techniques, the periodicity of the state can be extracted 
from the output state after applying the QFT, which exhibits the remarkable power of 
the QFT to determine the periodicity. Though the experimental imperfections cause the 
reduction of the fidelity, the powerful ability of quantum computers have been displayed. 
The QFT can be applied to many quantum algorithms and other quantum dynamics study. 
Recently, it has been used to implement the ordering-finding algorithm p71 and Shor's 



algorithm in NMR [28 
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Table 1. Measured NMR parameters for alanine dissolved in D 2 on a Bruker ARX500 



spectrometer. 



nuclei 


u/Hz J c '/Hz J c <*/Hz J c p/Hz Jh/Hz 


a 


-4320 


34.94 -1.2 5.5 


C a 


34.94 


53.81 143.21 


c p 


15793 -1.2 


53.81 5.1 


H 


1550 5.5 


143.21 5.5 



Table 2. The pulse sequence used to obtain the pseudo-pure state pooo = I?!^^ fr° m 
the thermal state, p eq = I lz + I 2z + h z . The refocusing -n pulses were omitted during the 



evolution th— due to Ji 3 << J 



(1) 
(2) 
(3) 
(4) 



2Jl2 



'12- 



Pulse Sequence for 1*1% 1% 



X 2J 12 
2 

X 4J23 



2Jl2 



- G, 



-G z 



Table 3 The pulse sequence for creating the labeled pseudo-pure state p oo — hz^i^2^3 
from the thermal state, p eq = Iq z + I± z + I 2z + hz- The refocusing n pulses were omitted. 



Pluse Sequence for Iq z I^ I 2 if 



(1) 
(2) 
(3) 
(4) 



7T 


1 


7T 


2 


7T 


3 


2_ 


y 


_2_ 


y 


_2_ 





G, 



- 


_j 


7T 


- s/ 


2 Joi 


4_ 


1° 


_J 


7T 


- y 


2 Jo2 


4_ 


1° 


_J 


7T 


- V 


2^03 


4 
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Figure Captions 



Fig.l Efficient circuit for three qubit quantum Fourier transform. H represents the 
Hadamard gate and conditionality on a spin being in the |1) state is represented by a filled 
circle on its time line. The final symbol inside the box represents the swap . 

Fig. 2 Experimental NMR spectra at different stages in the computation, (a) Experimen- 
tally measured thermal equilibrium spectra, acquired after a read-out pulse (|) on spin 
i. (b) Experimental spectra for the pseudo-pure state If . As desired, only one line 
is retained in each multiplet. (c) and (d) Output spectra of the QFT for the initial state 

= |0) + |2) + |4) + |6) and |0) + |4) , respectively. The positive (negative) lines correspond 
to a spin in |0) (|1)) before the readout pulse. The positive and negative lines counteract 
each other due to a spin in both |0) and |1) , so that no line is retained on the spin, (c) and 
(d) have been magnified by a factor of two for clarity. 

Fig. 3 Experimental deviation density matrices: (a) for the pseudo-pure state 1 000) , and 
for the states obtained after applying the QFT to the initial state \ip) of, (b) |0) + |2) + |4) + |6) , 
and (c) |0) + |4) , respectively. The rows are enumerated in the standard computational basis. 
The left and right columns denote the real and imaginary components (in arbitrary units), 
respectively. 

Fig. 4 Experimental NMR spectra for the labeling spin I at different stages during the 
computation. All spectra were recorded by a readout pulse n/2 on I . (a) Experimentally 
measured thermal equilibrium spectra, (b) Experimental spectra for the pseudo-pure state 
IqzI^I^I^- Only the transition Iq x Ii was presented, corresponding to the pure state 
Jf/^/g . (c) and (d) Output spectra of the simplified QFT for the initial state \ip) = |0) + 
|2) + |4) + |6) and |0) + |4) , respectively. The spectra of (b), (c) and (d) were added up 8 
times by phase cycling, (c) and (d) were magnified by a factor of two and four for clarity, 
respectively. 
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